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1. INTRODUCTION 

The HP-9825A calculator programs to be described here provide the 
capability to generate displays, such as those shown in Figures 1 and 2, 
which depict the dynamic pressure fields generated by Orbiter RCS thruster 
firings. Figure 1 was generated by the Jet Contour Plotting Program 
(JCPP), which is a stand-alone program that can be used to generate dynamic 
contours for an isolated RCS jet. The Orbiter Profile/Jet Contour Plotting 
Program (OP/JCPP) provides the capability to superimpose the plume contours 
for specific jets or jet clusters on front and side views of the Orbiter 
profile. The code for the JCPP and the OP/JCPP resides in Files 0 and 1, 
respectively, of the program tape. The tape is enclosed in a standard 
HP-9825A data cartridge. 

Both programs utilize the HP-9862A plotter and the HP-9866B thermal 
printer. The following four Read-Only Memory (ROM) blocks must be instal- 
led in the calculator before either program will load properly. 


ROM 

Part No. 

String-Advanced Programming 

9821 OA 

Matrix 

. 821 1A 

9885M Flexible Disk Drive 

982 17A 

9862A Plotter - General I/O - Extended I/O 

9821 4A 
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Figure 2. Typical OP/JCPP Output 
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2. PLUME MODELS 


Three plume modeling options are provided in the JCPP and the OP/UCPP. 
They are described separatf y in Sections 2.1 through 2.3. All of the 
models are axi symmetric and define dynamic pressure (q = 1/2 py ) as a 
function of a radial distance r and an angle e measured in a coordinate 
system which is illustrated in Figure 3. The X axis is coincident with 
the jet thrust line, and the origin of the system lies at the center of 
the nozzle exit plane. In the equations which follow, all distances are 
measured in feet, angles in degrees, velocities in feet per second, densi- 
ties in slugs per cubic foot, and pressure in pounds per square foot. 

2.1 MD/ECL MODEL 

This model (Reference 1) represents a curve fit of data published in 
Reference 2. Empirical equations for effluent density (p) and velocity (v) 
are defined separately for the continuum and the free-molecular flow 
regimes. In each flow regime, the flowfield is divided into two regions 
( | e j <_ 38.6° and |e| >_ 38.6°). As illustrated in Figure 4, the entire 
flowfield is thus divided into four areas, each having a particular set of 
density and velocity equations. These equations are shown in Table 1. 

To draw the contour line for a given value of q, arbitrary values are 
assigned to e and the appropriate equations in Table 1 are solved to find 
r. The solution is explicit in the free-molecular regime where v and W 
are functions of e only. In the continuum regime, v is a function of r» 
and an iterative solution is required. 

In the computational procedure, free-molecular flow is assumed 
initially, and the explicit solution 

r = (Wv 2 /2q) 1/2 

is calculated. The density is then calculated from 

p = W/r 2 . 

If p <_ 1.91885 x 10" 10 , the solution is complete. Otherwise, the value of 
r obtained from the free-molecular solution becomes the first guess in a 
Newton-Raphson iteration which uses the analytic derivative 
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Figure 4. MD/ECL Flow Regimes and Regions 
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for calculating corrections to r in the continuum regime. 

2.2 E&O MODEL 

In this model (Reference 3), the radial distance for given values of 
q and 0 is computed directly from the equation 

r = [546 f(e)/q] 1/2 


where 


f(e) = cos 7,14 (900/98.6) if e < 55.05 


and 


f(e) = 


0.04105 


exp [0.139934 (0 - 55.05)] 


if 0 > 55.05. 


The constants in the above equations result from assigning the following 
values to the parameters delineated in Reference 3: 


A = 1.1911 

Y = 1.28 

= 22032 
0 

r e ^ = 0.4 
A e /A* = 22. 

2.3 MD/QLDP MODEL 

This plume model is the one used in the MDAC Quick Look Disturbance 
Program (Reference 4). For given values of q and 0 , the radial distance is 
computed from the equation 

r = |[625 cos 8,47 (900/145)] /q j^. 

For this model, the magnitude of the angln 0 is limited to 145 degrees. 


8 


3. SETTING PHYSICAL PLOT BOUNDARIES 


The HP-9862A plotting board is a rectangle whose dimensions are 
approximately 17 inches by 12 inches. The plotter]? axis is parallel to 
the longer side and the 7 axis is parallel to the shorter side. The maxi- 
mum area accessible to the plotter's inking pen is a smaller rectangle whose 
dimensions are 15 inches by 10 inches. A plot can be drawn on any sheet of 
paper small enough to fit on the plotting board. T he pape r is held in place 
by electrostatic force that is controller oy the QHART H0U}) button on the 
plotter control panel. As indicated in Figure 5, a small sheet of paper 
can be located at an arbitrary position on the board, and the plotting area 
can be bounded on the paper by a simple procedure that is described in the 
next paragraph. 


For convenience, fiducial marks should be inscribed on the sheet of 
'lotting paper at the lower left and upper right corners of the desired 
plotting area before it is placed o n the {j otting board. After the paper 
h as been secured by pres sing th e (£HART HOLp) button, the user presses the 
J»EN~ Up land then the (LOWER LEF^) buttons on the plotter control panel. 
After the pen comes to rest, its position is adjusted to coincide with the 
lower left cor ner of the desired plotting area by rotating two knobs adja- 
cent to the (L OWER LEF T) button. These two knobs move the pen parallel to 
the plotter 7 and Taxes independently. After the pen point has been posi- 
tioned over the lower left fiducial mark as closely as the eye can judge, 
the (P EN DOW fj) and ( Ten UP ^ buttons can be pressed in sequence to mark the 
exact pen position on the paper. If necessary, fine adjustments of the pen 
position can then be made. After the user has adjusted the pen position 
at the lower left corner to his satisfaction, he presses the CUPPER RIGHT 
button and positions the pen over the upper right fiducial mark by using 
the adjacent control knobs. Again, the exact position of the pen can be 

buttons ir. succession 


determined by pressing the (^P EN DOW ?!) and( j>EN UP, 
to produce a dot of ink on the paper. 


The sequence of operations described in the preceding paragraph is 
important. The adjustment of the lower left corner must always be made 
first. If for any reason a readjustment is made at the lower left corner, 
it must be iorowed by a readjustment at the upper right corner. 
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PEN'S LIMIT OF ACCESSIBILITY 
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Figure 5. Physical Plot Boundaries 





4. JCPP PROGRAM 


4.1 OPERATING INSTRUCTIONS 

1. Turn on power switches on the HP-9825A calculator, the HP-9862 
plotter, and the HP-9866B printer. 

2. Insert data cartridge containing program tape into the slot to 
the left of the display window on the HP-9825A calculator. 

3. Pres s( j*E W I NjT ) * . 

4. Type in the instruction "ldm0" + (load memory from tape file 0). 
Verify that the instruction appears correctly in the calculator 
display window, then press (E XECUT E) . This will cause the pro- 
gram instructions and the standard input data base to be loaded 
into memory. 

5. When th e lazy T ("i-") character appears in the display window, 
press CREW I N^ ) . The data cartridge may be removed when the 
amber tape activity light below the cartridge slot ceases to glow. 

6. Place a sheet of paper on the HP-9862A plotting board and adjust 
the physical bounds of the plotting area as explained in Section 3. 

7. Press (^ CMTI NUE) . The calculator will respond by listing the 
current contents of the input data base (see Section 4.2) on the 
HP-9866B printer. 

8. The query "WHAT REGISTER?" should now appear in the display 
window. If no ( further ) changes to the input data base are 
desired, press CcONTINui) and skip ahead to Step 11. Otherwise, 
go to Step 9. 


* Individual keys on the calculator keyboard a re des ignated by circumscribing 
the characters appearing on the key (e.g. , (^ EXECUTE ) ). 

J* 

Instructions or messages appearing in the calculator display window are 
enclosed in quotation marks. 
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9. Enter the number of the data register (whos e conten ts are to be 
changed) into the display, and then press C ONTINU E) . The cal- 
culator will respond by printing the current contents of the 
specified register* on the HP-9866B printer. 

10. The query "NEW VALUE?" should now appear in the dis 
Enter the desired value + into the display and press 
The calculator will print the new value on the HP-9866B. Go back 
to Step 8. 

11. The program will now begin automatic execution and will continue 
plotting until all of the specified dynamic pressure contours 
have been drawn. If any input parameters were changed, the 
revised input data base will be listed in its entirety before 
plotting begins. 

12. When plotting is complete, the calculator will issue a beep sound 
and display the message "FINIS" in the display window. If a new 
set of contours is to be drawn on the same sheet of paper, go back 
to Step 7. Otherwise go to Step 13. 

13. Press the (C HART HOL D) button on the plotter control panel to 
release the electrostatic force and then remove the completed 
chart. If a new chart is to be prepared, go back to Step 6. 
Otherwise, go to Step 14. 

14. Turn off the power switches on the HP-3825A calculator, the HP-9862A 
plotter, and the HP-9866B printer. 


pl ay windo w. 
CONTINUE. 


* If the operator's entry was not a valid input data register number, nothing 
will be printed. In such a case, go back to Step 8. 

*t* 

The entry may be a number or any valid expression. If an expression is 
entered, the calculator will evaluate it and store the result in the speci- 
fied input data register. If no entry is made, the contents of the input 
data register will not be changed. 
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4.2 INPUT DATA 


The standard input data base which resides on the program tape is 
shown in Table 2. The operator may change any or all of the input values 
at the beginning of a run by use of the editing logic described in Section 
4.1. If another run is initiated without reloading memory from the program 
tape, the previously-edited data base becomes the starting point for the 
new run. Individual input parameters are explained in Sections 4.2.1 
through 4.2.5. 


4.2.1 

Logical Plot Boundaries 

(Registers 

0 - 3) 

FT 

-6. 0 0 0 0 0 0 0 0 E 01 

0 : ';■■■ 

M I H 

FT 

5, 09000000E 02 

l: 1 

MAX 

FT 

- 1 . 8 0 0 0 0 0 0 0 E 0 2 

• 'i 

MIN 

FT 

1 . 80000 000 E 02 

■~i . 

o » i 

MAX 


The logical boundaries of the plotting area are defined in terms of 
the Cartesian flowfield coordinates X = r cos e and Y = r sin e (see 
Figure 3). The logical coordinates (X mi - n , Y^j correspond to the lower 
left corner of the physical boundary described in Section 3. The logical 
coordinates /X , Y \ correspond to the upper right corner of the 
physical boundary. 

The logical and physical boundaries can be adjusted to plot any 
portion of the flowfield at any desired scale. However, the equation 



AX. X max " X min 


must be satisfied if distortions are to be avoided, where aY and aT are 
as defined in Figure 5. The standard input values for registers 0-3 
yield a scale of 1" = 40* when a)T = 14" and aY = 9". 

4.2.2 Tic Interval (Register 4) 

FT 2 . 0009000UE 0i 4: TIC INTERVAL 

If register 4 contains a number greater than zero, the flowfield 
coordinate axes will be drawn on the display, with tic marks at the speci- 
fied interval. If register 4 contains a zero or a negative number, the 
coordinate axes will not he drawn. 
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Table 2. JCPP Standard Input Data Base 


LU 

I! CO 


d* LU 

LU 
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Register 4 should not be assigned a positive value if the logical 
plot boundaries (Section 4.2.1) are defined such that the origin of the 
flowfield coordinate system lies outside the plotting area. If a positive 
value Is assigned in such a case, an error condition will result when the 
plotter attempts to draw the axes. 

4.2.3 Plume Model (Register 5) 

0 5 : P L. U M E N 0 D E L. •" u - N D . " E C L. < 1 - E D » 2 = M D 

Any of the plume models described in Section 2 can be utilized by 
loading the appropriate value into register 5. 

4.2.4 Contour Specifications (Registers 6-8) 

L B / F T 1 2 1 . 0 u 0 6 0 0 0 0 E - 0 5 € : Q M I N 

i . 0001000 00 E 01 7: Q RflTf'O BETWEEN CONTOURC 

6 8: HO. CONTOURS 

The value assigned to register 6 defines the dynamic pressure for the 
outermost flowfield contour. For each succeeding contour, the dynamic 
pressure is multiplied by the contents of register 7. Register 8 defines 
the total number of contour lines to be drawn. 

4.2.5 Jet Multiplier (Register 9) 

1 9 ‘ NO. JETS 

Provisions are made for approximating the flowfield of clustered jets 
firing simultaneously. When this option is used, the number of jets is 
entered in register 9. It is assumed that the dynamic pressure in an 
N-jet flowfield is N times greater at every point than the dynamic pressure 
produced by a single jet . No representations are made as to the validity 
of the assumption; certainly it is not accurate at short distances from 
the jet exit(s). It is left to the user to determine how much credence 
should be given to multi-jet data generated in this manner. 
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5. OP/JCPP PROGRAM 


5.1 OPERATING INSTRUCTIONS 

1. Turn on power switches on the HP-9825A calculator, the HP-9862A 
plotter, and the HP-9866B printer. 

2. Insert data cartridge containing program tape into the slot to 
the left of the display window on the HP-9825A calculator. 

3. Press C jEWIND ^*. 

4. Type in the instruction "ldml" f (load memory from tape file 1). 
Verify that the instruction appears correctly in the calculator 
display window, then press ^EXECUTE) . This will cause the pro- 
gram instructions, program constants, and the standard input data 
base to be loaded into memory. 

5. When th e lazy T ("i-") character appears in the display window, 
press( ^REWIND ^. The data cartridge may be removed when the amber 
tape activity light below the cartridge slot ceases to glow. 


6. Place a sheet of paper on the HP-9862A plotting board and adjust 
the physical bounds of the plotting area as explained in Section 3. 

7. Press (^C ONTINUE ^) . The calculator will respond by listing the 
current contents of the input data base (see Section 5.2) on the 
HP-9866B printer. 

8. The query "WHAT REGISTER?" should now appear in the display 
window. If no (f urther ) changes to the input data base are 
desired, press C ONTINU E]) and skip ahead to Step 11. Otherwise, 
go to Step 9. 


* Individual keys on the calculator keyboard a re desi gnated by circumscribing 
the characters appearing on the key (e.g. , ( JXECUTE ^ ) . 

Instructions or messages appearing in the calculator display window are 
enclosed in quotation marks. 
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9. Enter the number of the data register (whos e conten ts are to be 
changed) into the display, and then press ( ^NTINUE ^ . The cal- 
culator will respond by printing the current contents of the speci- 
fied register* on the HP-9866B printer. 


10. The query "NEW VALUE?" should now appear in the displa y windo w. 
Enter the desired value + into the display and press (^ CONTINU E) . 
The calculator will print the new value on the HP-9866B, Go back 
to Step 8. 

11. If any input parameters have been changed, the program will list 
the revised input data base in its entirety. 

12. The query "0V PROFILE?" should now appear in the display window. 

If it is not desired to have the plotter draw the profile view(s) 
of the Orb ite r vehic le, load a zero ("0") i nto the display and 
then press (^ CONTINU E]) . Otherwise, press (^C ONTINU E]) without 
making any entry, or else lo ad any non-zero value into the dis- 
play and then press ( ^ONTINU ?) . 

13. The program will now begin automatic execution and will continue 
plotting until all of the specified dynamic pressure contours have 
been drawn. If no entry was made in Step 12, or if a non-zero 
value was entered, the profile view(s) of the Orbiter will be 
drawn before the dynamic pressure contours are plotted. 


14. When plotting is complete, the calculator will issue a beep sound 
and display the message "FINIS" in the display window. If a new 
set of contours is to be drawn on the same sheet of paper, go 
back to Step 7. Otherwise, go to Step 15. 


* If the operator's entry was not a valid input data register number, nothing 
will be printed. In such a case, go back to Step 8. 

+ The entry may be a number or any valid expression. If an expression is 
entered, the calculator will evaluate it and store the result in the speci- 
fied input data register. If no entry is made, the contents of the 
input data register will not be changed. 


17 


15, Press the (C HART HOL D) button on the plotter control panel to 
release the electrostatic force and remove the completed chart. 

If a new chart is to be prepared, go back to Step 6. Otherwise, 
go to Step 16. 

16. Turn off the power switches on the HP-9825A calculator, the HP-9862A 
plotter, and the HP-9866B printer. 

5.2 INPUT DATA 

The standard input data base which resides on the program tape is 
shown in Table 3. The operator may change any or all of the input values 
at the beginning of a run by use of the editing logic described in Section 
5.1. If another run is initiated without reloading memory from the program 
tape, the previously-edited data base becomes the starting point for the 
new run. Individual input parameters are explained in Sections 5.2.1 
through 5.2.5. 

5.2.1 Logical Plot Boundaries (Registers 0-5) 

FT 
FT 
FT 
FT 
FT 
FT 

The logical plotting boundaries are defined in terms of Orbiter body 
axes (OB) coordinates. This coordinate system is described in Figure 12 
of Reference 5. 

Figure 6 shows the relationship between the logical plot boundaries 
and the physical boundaries established by the procedure explained in 
Section 3. As indicated, the plotting area is partitioned to provide a 
side view and a front view of the Orbiter. The location of the partition 
and the scale of the plotted figures can be adjusted by changing the values 
in registers 0-5. The front view can be suppressed by setting YQ Bmflx = 
Y 0Bmin* or the s1de view can be suppressed by setting X 0Bmax = X 0Bmi|) . In 
any event, the equation 

a 7 PoBmax " ^OBminj 

^ Y 0Bmax ~ Y 0Bmin^ + ^ X 0Bmax ” *0Bmin^ 


1 . 2000000 0E 02 0 

YOB Mill 

2 „ 00 000 000 E 02 1 

2 ob nh: . 

1 -2000000 0E 02 2 

YOB Mil! 

1.20000000E 02 3 

yob hr:-: 

2. 8000000QE 02 4 

ZOB MIN 

O.00000000E 01 5 

ZOB MAX 
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Table 3. OP/JCPP Standard Input Data Base 
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Figure 6. OP/JCPP Logical Plot Boundaries 



must *'** satisfied if distortions are to be avoided. The standard input 
values for registers 0-5 yield a scale of 1" « 40' when a 7 = 14" and 
aT = 9". 

3.2.2 Orbiter CG Location (Registers 6-8) 

sJn 1.0767000OE S3 6: OV CG STR 

IN 0.00Q00000E 00 7: OV CG UL 

IN 3.74100000E 02 8: OV CG Wl 

The Orbiter center of gravity (CG) location, which is the origin of 
the OB coordinate system, is defined in terms of Orbiter structural body 
coordinates. The structural body coordinate system is defined in Figure 23 
of Reference 5. In OP/JCPP, the X, Y, and Z structural body coordinates 
are designated, respectively, as the station (STA), buttock line (BL), 
and water line (WL). 

5.2.3 Plume Model (Register 9) 

0 9s PLUME MODEL 1 0-MB' ECL* 1=ELD» 2=MD / ’QL 

Any of the plume models described in Section 2 can be utilized by 
loading the appropriate value into register 9. 

5.2.4 Contour Specifications (Registers 10 - 12) 

LB--FT12 2. 50O000O0E-02 10: Q MIN 

2.00000000E 00 11 : Q RATIO BETWEEN CONTOURS 

5 12: NO. CONTOURS 

The value assigned to register 10 defines the dynamic pressure for the 
outermost flowfield contour. For each succeeding contour, the dynamic pres- 
sure is multiplied by the contents of register 11. Register 12 defines the 
total number of contour lines to be drawn for each active jet or jet cluster. 


5.2.5 Jet Multipliers (Registers 13 - 22) 
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JETS 
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JETS 
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NO. 
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LH 

♦2 JETS 

' 09 1 f 2 > 

-- 
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NO. 
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RH 

♦ IETS 

( 0 » 1 * 2 * 3 ) 
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Registers 13-22 contain the number of active jets in each of 10 RCS 
jet clusters. Each cluster is identified in terms of its location (e.g. t 
AFT LH) and the general direction of the force it exerts on the Orbiter 
(e.g., +X), as measured in the Orbiter body axes coordinate system. The 
direction of the jet effluent is opposite to the direction of the force it 
produces on the Orbiter. Downward-firing (-Z) jets are not modeled in the 
OP/JCPP. 

A value of zero in one of the jet multiplier registers signifies that 
all the jets in the appropriate cluster are inactive, and that no contours 
are to be drawn for that cluster. For the purpose of generating the con- 
tours for any particular cluster, it is assumed that the dynamic pressure 
in an N-jet flowfield is N times greater at every point than the dynamic 
pressure produced by a single jet . The possibility of interactions between 
the flowfields produced by different jet clusters is ignored. The effect 
of Orbiter structure on the flowfields is also ignored. No representations 
are made as to the validity of these assumptions and simplifications; it 
is left to the user to determine how much credence should be given to the 
flowfield data thus generated. 

The origin of each flowfielr is taken to be at the mean location of 
the centers of the exit planes for the i 'i vi dual jets in the appropriate 
cluster. The origin and orientation of each fowfield is defined in Table 4. 
The information given there is based on data outained from Reference 6. 
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Table 4. RCS Jet Cluster Flowfleld Data 
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